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Last year, we reported our first results on the determination of Gasser-Leutwyler coefficients using partially 
quenched lattice QCD with three flavors of dynamical staggered quarks. We give an update on our progress in 
determining two of these coefficients, including an exhaustive effort to estimate all sources of systematic error. At 
this conference, we have heard about algorithmic techniques to reduce staggered flavor symmetry breaking and a 
method to incorporate staggered flavor breaking into the partially quenched chiral Lagrangian. We comment on 
our plans to integrate these developments into our ongoing program. 



1. Introduction 

The broad application of effective field theories 
like the chiral Lagrangian to the description of 
low energy phenomena in QCD has been one of 
the major theoretical advances of the last quar- 
ter century pj] . One substantial limitation of this 
approach is that the chiral Lagrangian has addi- 
tional symmetry H that precludes the determi- 
nation of all the low energy constants (LEC's) of 
the effective theory using only experimental in- 
put. This problem is resolved by computing low 
energy observables directly in QCD and fixing the 
undetermined LEC's by matching to predictions 
of chiral perturbation theory (ChPT). 

Unfortunately, the low energy regime of QCD 
is non-perturbative, so relevant observables are 
not so easily computed. This has led to an entire 
industry of estimating LEC's (with uncontrolled 
systematic errors) using model calculations and 
phenomenology, culminating in a generally ac- 
cepted set of values [||. Lattice QCD has also 
made some inroads in this area, with its own set of 
systematic errors, including errors due to Nf =/= 3 
and dynamical quarks that are too heavy. See 
Wittig's review talk for a current summary Q . 

Recently it was realized that ChPT and the 
corresponding partially quenched chiral pertur- 
bation theory (pqChPT) [§,|| share the same low 
energy constants As partially quenched 

QCD (pqQCD) simulations are somewhat less 
computationally demanding, recent efforts have 
been made to determine the LEC's, particularly 



those affecting the mass of the up quark || . How- 
ever, as the number of light dynamical sea quarks 
Nf doesn't appear explicitly in the chiral La- 
grangian, LEC's computed with Nf ^ 3 have no a 
priori relation to those of phenomenological inter- 
est, although a posteriori evidence suggests that 
the Nf dependence may be mild. Here, we will 
review our recent Nf = 3 calculations JTo) , em- 
phasizing our program to estimate all expected 
systematic errors, and discuss our plans for re- 
ducing the systematic errors in the future. 



2. LEC's and systematic error 

Our goal is to determine the low energy con- 
stants 2a§ — as and 2a^ — a± [] by studying the 
sea and valence quark mass (mj and my) depen- 
dence of the pseudo-Goldstone meson mass M v 
and decay constant NLO pqChPT predicts 
the dependence in terms of the LEC's, with loops 
cut off at A x = 47r/, as 
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2/i/A?. is related to the chiral con- 



densate fj, and M refers to a QCD meson mass, 
without QED contributions. Similar formulae for 
fixed ms can be used to extract only a§ and 
2a% — a 5- We refer readers first to our earlier 
work [|lOj and focus on details not covered there. 

The primary sources of systematic error are 
finite volume, finite lattice spacing, finite step 
size, flavor (or "taste" |nj) symmetry violations 
and the choice of sea and valence quark masses 
for observables included in the fitting procedure. 
The errors are not independent. For example, as 
smaller quark masses are included in the analysis, 
larger volumes are required to minimize finite vol- 
ume errors. Also, at smaller lattice spacings all 
systematic errors due to lattice artifacts are re- 
duced, including flavor symmetry violations. So, 
while not justified by the preceding observation, 
we have chosen to make generous estimates of our 
systematic errors and add them in quadrature. 

Our procedure for estimating systematic errors 
is simple. Holding all other parameters fixed, we 
independently vary the volume, step size, lattice 
spacing and sea quark mass. To estimate fitting 
systematics, we vary the range of quark masses 
included in the fit, ensuring that all masses are 
within the expected radius of convergence of the 
effective theory. To estimate flavor symmetry 
breaking still remaining at a given lattice spacing, 
we perform HYP blocking |l2| of the configura- 
tions and recompute the observables. This pro- 
cedure may seem ad hoc but it certainly has the 
correct continuum limit as the blocking becomes 
irrelevant for sufficiently smooth gauge fields. 

As our paper shows, this recipe for estimat- 
ing systematic errors works well with M T . HYP 
blocking doesn't significantly change the leading 
order behavior, as expected since local blocking 
shouldn't effect quantities measured at long dis- 
tances on the lattice. In the presence of flavor 
breaking, the pions contributing to the chiral logs 
are too heavy and systematically shift NLO terms 
away from their continuum values. With HYP 



blocking, these pions become more degenerate 
and the systematic error is reduced. 

An interesting contrast is what happens to f v . 
Since the dominant contributions are local, Fig. 
clearly shows the effect of HYP blocking on the 
leading order behavior of f n . This seems dis- 
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Figure 1. f n with and without HYP blocking on 
16 3 x 32, 13 = 5.3, and m s = 0.01. 

concerting at first glance until one notices the 
unimproved staggered /„. is substantially higher 
than the PDG value whereas the HYP blocked 
/„■ is consistent with it. The real difficulty comes 
when extracting 015: for unimproved staggered 
a 5 = 2.134(17) and for HYP a 5 = 0.242(34). 
It doesn't seem reasonable to treat the differ- 
ence in the central values as an estimate of sys- 
tematic error. Luckily, a§ is not affected by the 
Kaplan-Manohar ambiguity so phenomenological 
estimates are likely to be more reliable. On the 
other hand, these estimates currently have 20% 
errors so we would like to consider how to better 
determine LEC's like as in the next section. 

3. Future plans 

Our original calculation jnj centered around a 
simulation of a single dynamical sea quark mass 
at a" 1 « 1.3 GeV, thus it was not possible to 
extract 2a§ — 0:4. However, our study included 
supporting simulations with several different sea 
quark masses at a -1 w 0.7 GeV on comparable 



physical volumes. So, we offer Fig. | as a demon- 
stration of the technique to extract 2a§ — and 
2a% — a$ simultaneously. This figure shows a sin- 
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Figure 2. i? M = M%(m s )m v / M%(mv)m s on 
8 3 x 32 lattices at a" 1 w 0.7 GeV 

gle fit to four simulations with x 2 /dof ~ 1 and 
should be compared to the corresponding figure 
in Jl0| where each simulation is fit independently. 
We expect to perform a detailed extraction of 
2cts — cti by this method, including a study of 
systematic errors, when sufficient data becomes 
available at smaller lattice spacings. 

In section |], we argued that HYP blocking is 
useful for estimating systematic errors due to fla- 
vor symmetry breaking. At this conference, we 
have heard about recent efforts to find an efficient 
dynamical HYP algorithm |H|[l4j. We hope to 
generate an ensemble of Nf — 3 dynamical HYP 
lattices which will be used to calculate kaon weak 
matrix elements for e'/e jl5| as well as Gasser- 
Leutwyler coefficients. This should help resolve 
such questions as how to estimate the systemat- 
ics of flavor breaking in quantities like / T . 

Also at this conference, Bernard provided an 
update on his program to incorporate flavor sym- 
metry breaking into ChPT pj| . This is an im- 
portant development for our program. It alto- 
gether eliminates flavor breaking as a source of 
systematic error as the effects are correctly mod- 
eled in our fitting function. It may also increase 
the range over which the data are well represented 



by our fitting functions (still within the expected 
radius of convergence, of course) . 

Finally, we note that while our determination 
of 2a% — a§ is completely dominated by system- 
atic error, its systematic error is not the dominant 
source of uncertainty in the ratio m u /md- The 
dominant uncertainty comes from the unknown 
electromagnetic contributions to the masses of 
the pseudoscalar mesons. A technique was de- 
vised to measure these effects in quenched lattice 
simulations []l6| and was presented at Lattice '96 
pTy . We intend to explore the extension of this 
work to simulations with dynamical fermions and 
hope to present our first results at Lattice '03. 
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